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Abstract 

We apply a new anticommuting path integral technique to clarify the fermionic struc- 
ture of the 2D Ising model with quenched site dilution. In the iV-replica scheme, the 
model is explicitly reformulated as a theory of interacting fermions on a lattice. An 
unusual feature is that the leading term of interaction in the exact lattice theory is 
of order 2N in fermions, where N is the number of replicas. The continuum-limit 
approximation near T c for weak dilution produces, however, an effective four-fermion 
interaction. In particular, this implies the doubled-logarithmic singularity in the spe- 
cific heat near T c for weak site dilution. The exact value of the initial slope is also 
obtained. 
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1 Introduction 



From the experimental point of view, the site dilution provides, probably, the most simple 
way to realize quenched disorder in real magnetic materials. In this case, some amount of 
the magnetic atoms in a sample, chosen at random, are replaced by nonmagnetic impurity 
atoms. Another sort of disorder is bond dilution, merely preferred in theoretical studies, in 
which case some of the lattice bonds are assumed to be broken. The two-dimensional Ising 
model (2DIM) is a natural object to analyze the effects of disorder since it admits the exact 
analytic solution in the pure case Q. This model, on the other hand, is too delicate to make 
definite predictions on common phenomenological grounds. In particular, from the point of 
view of the Harris criterion [fj, |3| (which states that the critical behaviour is modified by 
quenched disorder if a > 0, and is unchanged if a < 0, where a is the specific heat exponent 
for the pure system) the 2D Ising model is marginal, since a = in the pure case. The 2D 
Ising model with quenched dilution has been intensively analyzed during the last decades 
both by theoretical methods PI-fllBf and in the precise Monte-Carlo experiments [l7|-[22 



The theoretical studies here were merely concerned with bond dilution, making use of the 
fermionic interpretation of the Ising model ]Pfl-||16||. The case of site dilution, however, has 
not yet been properly analyzed within the fermionic approach. The site dilution can be 
viewed as locally correlated bond dilution. This makes the traditional transfer-matrix and 
combinatorial methods of the fermionization, which is the basis for theoretical analysis, less 
suitable in the case of site disorder Meanwhile, the Monte-Carlo studies [|20| and related 



analyses |23] give evidence for a different behaviour, at least for moderated and strong site 
dilution, from that predicted in the theoretical studies of weak bond dilution || 14j. The 
alternative theories for bond dilution are also to be taken into account jTT[ |T2[. A number 
of questions thus arise about the universality between site and bond dilution, and between 
weak and strong dilution in both the cases ||22||-||25||. Also, see for a critical review. 

In this Letter we apply a novel anticommuting path integral technique to clarify the 
fermionic structure of the two-dimensional Ising model with quenched site dilution. At the 
first stage, we transform the partition function of the 2D Ising model with fixed site dilution 
into a Gaussian fermionic integral. The averaging over the disorder in the iV-replica scheme 
then yields a lattice theory with compact multifermion interaction. An interesting feature 
is that the term of interaction here is of order 2N in fermions, where N is the number 
of replicas. Notice that the formal limit N — * is assumed at final stages. For weak 
dilution, however, the effective continuum- limit theory responsible for critical behaviour 
near T c appears to be the standard N = Gross-Neveu model with four-fermion interaction, 
already well known in the DD-SSL theories for weak bond dilution || [14]]. In particular, 
this implies the doubled-logarithmic singularity in the specific heat and log-corrections in 
other thermodynamic functions near T c for weak site dilution. An effective four-fermion 
term arises in the continuum theory, in our case, due to the interplay of the short-wave and 
long-wave lattice fermionic modes. The parameters of the N = Gross-Neveu model for site 
dilution are then explicitly evaluated. These parameters include some characteristic lattice 
fermionic averages for the pure Ising model at T c . The exact value for the initial slope of 
the T c (p)-p curve is also obtained. The basic results are given in (H), (|T0|) , (|15|), fl3l[), ([37]). 
Further comments are added in the following sections. 



1 



2 The model 



Let us start with the two-dimensional Ising model with fixed distribution of the nonmagnetic 
sites over the lattice (fixed site dilution). At the first stages, we also assume arbitrary 
inhomogeneous distribution of the coupling parameters over the lattice bonds since this will 
not be an obstacle for fermionization. The hamiltonian is: 

/3H X/ [ bm+ln ymnym+ln&mn&m+ln ymnymn+l&mn&mn+l ] j (1) 

ran 

where the Ising spins, <J m n — 3re disposed at the lattice sites, twi, with 771, n — 1, 2, L , 
L is the lattice length. To introduce site dilution, we supply each Ising spin by random 
variable y mn = 0, 1 playing the role of the effective magnetic moment of the Ising spin. If 
y mn = 1, the given site is normal, if y mn = 0, the site is dilute. Finally, = /3J^°$ are 
the dimensionless bond couplings, are the exchange energies, (3 = 1/kT is the inverse 
temperature in the energy units. > corresponds to the ferromagnetic bond. For fixed 
disorder, the partition function and free energy are: 

Z{y} ^e-W'I'^e-W, (2) 
to 

where the sum is taken over all the possible spin configurations provided by a mn = ±1 at 
each site. Taking into account the identity for a typical Boltzmann weight: exp {byy'aa') = 
cosh(&yy')(l + yy'aa' tanhfr), which readily follows from aa' = ±1, yy' = 0, 1, the partition 
function can be written in the form: Z {y} = R{y} Q {y}, where Q{y} is the reduced 
partition function and R{y} is a nonsingular spin-independent prefactor to be ignored in 
what follows. The reduced partition function is: 

Q {y} = Sp Y\( 1 + tln+ln ymnym+ln°'mn 'm+ln) (1 + ymnymn+l cr mn cr mn+l) j (3) 

to mn 

where = tanh and we assume a properly normalized spin averaging such that 
Sp (1) = 1, Sp (<r mn ) = at each site. 

Since we are interested in quenched disorder, the averaging over the random-site (RS) 
impurities is to be performed for the free energy rather than for the partition function itself. 
The standard device to avoid the averaging of the logarithm is the replica trick: 



\nQ{y}= Hm-(Q"M-1). (4) 

In this scheme, we take iV identical copies of the original partition function and average 
Q N {y}- The formal limit iV — > is assumed at final stages. We will also assume the 
following distribution of the RS impurities in the definition of the averaging: 

W(y mn ) =p6(l -y mn ) + (1 -p)S(y mn ) , (5) 

where S( ) are the correspondent Kronecker's symbols, p is the probability that any given 
site, chosen at random, is occupied by the normal Ising spin, while 1 — p is the probabilty 
that the given site is dilute. 
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3 Fermionization 



Before averaging over the disorder, we have to transform Q{y} into a fermionic Gaussian 
integral. Let us remember that Grassmann variables (nonquantum fermionic fields) are the 
purely anticommuting fermionic symbols. Given a set of Grassmann variables oi, a2, ■■■ , ajv, 
we have a^aj + = 0, a| = 0. The notion of the integral over Grassmann variables 
(fermionic path integral) was first introduced by Beresin [26|. The elementary rules of 
integration for one variable are [|26|: 



J ddj • dj = 1 , J ddj ■ 1 



0. 



(6) 



In a multidimensional integral, the differential symbols da\, da2, . . . ,da^ are again anticom- 
muting with each other and with the variables. The fermionic exponential is assumed in 
the sense of its series expansion. For a finite number of the variables, the exponential se- 
ries definitely terminates at some stage due to the property a 2 = 0. In the field-theoretical 
language, the fermionic form in the exponential is called action. 

The transformation of Q{y} into a fermionic integral can be realized by following literally 

For a 



the fermionization procedure first developed for the pure 2D Ising model p7f-p9[ 
recent comment see also |y], HU- The first step is the fermionic factorization of the local 
bond Boltzmann weights |27|, pjjfl . For the whole lattice, we introduce a set of the totally 
anticommuting Grassmann variables, a 



»>n, <d, b mn , b mn , and write: 



1 ~l~ ^m+lnZ/mn?/m+ln Cr »Tm cr r?i+ln 

"mn c V x ' "'rani) 



rln* An e , a mnO, mn 



)(l+t (1) 



m+ln ^mnVra+XrP ~m+ln) ' 



1 ~l~ ^mn+iyranDmn+l'^irin'^mn+l 

1 r> i _|_ Jj 



dh* dh p"" 1 """" 
uu mn LLU mn " 



) (1 + tlJ n+1 b* mn y 



mn+l L 'mnii 'mn+l" mn+l 



(7) 



These identities can be readily checked using the elementary rules like and taking into 
account that e aa * = 1 + aa*, since (aa*) 2 = 0. Neglecting the sign of the Gaussian fermionic 
averaging, the weights are now presented as y4 mn y4^ +ln , _B mr jS,^ n+1 , where the separable 
factors (to be called Grassmann factors) are to be identified from (^). At the next stage we 
group together, over the whole lattice, the four factors with the same Ising spin (the same 
index mn) and then sum over a mn = ±1 in each group of factors independently, thus passing 
to a purely fermionic expression for Q{y}. The above four factors with the same a mn come 
by factorization of the four bonds attached to the given mn site, they all include also the 
same variable y mn . There is also the ordering problem, to be solved, for the noncommuting 
Grassmann factors in their global products. This problem is related to the requirement 
that we can actually keep nearby the four factors with the same spin when averaging over 
a mn = ±1 at each site. In the 2D case, the ordering problem is resolved in terms of the 



'mirror-ordered factorization procedure' for the global products of weights [27]-[29 . The 



appearance (as compared with the pure case) of additional parameters y mn in the factors 
arising in (|?|) does not matter in this respect. The averaging over a mn ± 1 in the local group 
of four factors yields an even fermionic polynomial. This polynomial can be represented as 
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an exponential factor associated with the mn lattice site. This is just the exponential factor 
corresponding to the y mn term in (§D below. Thus we come to the result: 

Q {y} = / II db mn db mn da mn da mn exp ^ { a mn a mn + b mn b mn 

mn mn 

_l 7/ 2 \ n h -i-fWfWn* h* + rt* + / (2) fr* Ho +fr "ill fS^i 

1 smti L mn mn 1 "mti"mii u 'm-ln"mn-l 1 rran^m-ln 1 v mn mn—1) V^mn 1 u mnJ J J " \ u / 



In the RS pure case {y mn = 1) we come back to Eq. (11) of Ref. [£7] (inhomogeneous 
distribution of bond couplings was already assumed in p7| ). In fact, noting (^), one can 
guess (H) directly from Eq. (11) of Ref. p7j without special calculation. The partition 



function with fixed disorder is represented now as a Gaussian fermionic integral. 

The integral (§) can be in turn simplified integrating out the a mn , b mn fields by means 
of the identity / dbda exp(Aa& + aL + L b) = A + LL' = A exp (A^ 1 LL' ), where a, b 
are Grassmann variables, L, L are linear fermionic forms independent of a, b, and A is a 
parameter p2| . For A = the result is LL'. Integrating out the a mn ,b mn fields in this way, 



we obtain a reduced integral in terms of variables a mn , b mn . Changing the notation for the 
fields, for the sake of visual comfort, as follows: a mn , b* mn — > c mn , — c mn , we find the reduced 
integral in the form: 

0\v\= fV\dr dr v 2 exnV [w" 2 r r - ?/ 2 f (1) f (2) 

% l!»J / J. 1_ ^^mn^^mn i> mn H / y \i>mn^mn^mn i)mn mn v mr 



mn °mn°mn ^m-ln^mn-l 

mn 

_L ( r UfW, _ +(2) 7; \1 (q) 

1 V^mn 1 ^mnJ rmn ^m-ln mn L 'mn— 11 ' V / 



where y£ n exp (y m 2 n c mn c mn ) = y mn + c mn c mn , there is no true singularity at y 2 mn = 0. The 
fermionic integrals t^) and (|9|) are equivalent to each other and to (0). This equivalence 
holds even independently of the formal interpretation of the variables y mn . However, we have 
essentially used the property y mn = 0, 1 by passing from Z{y} to Q{y}. If y mn = 0, 1, then 
one can put y mn instead of y mn in the above integrals. In what follows we assume definitely 
that y mn = 0, 1. The structure of the integral can be illuminated then most obviously by 
expanding the local factor in (|9|) over the 'eigenstates' with y mn = and y mn = 1: 

Q{y}= Ii d c mn dc mn J] [5(1- y mn ) exp 

(S m n) + 5 (y mn ) c m n C mn > 

(10) 



■mn 



where 5(1 — y mn ) = y mn and 5 (ymn) = 1 — y mn are the correspondent Kronecker symbols, 
y mn = 0, 1, and S mn is the local density of action for the RS pure case, see fliTI). 

In the RS pure case, with y mn = 1 at all sites, we obtain: 



Q W = / II D mn exp ^2 Smn = / II dc mn dc mn exp ^2 [ 

mn mn mn mn 

+ ( c +c Ht (1) c 1 -t (2) c il-t (1) t (2) c ic 1 

~ y-*mn ~ '-"mn ) ^mn ™-!" v mn L "n^n—l J raii l 'mn t ™-l" l ™«-l 



where D mn and SVrm are the abbreviations for the local terms in the fermionic measure and 
the action, respectively. The inhomogeneous distribution of bond couplings is still preserved. 
Therefore, integral (|TT|) can be used as a starting point to analyze bond dilution. We can 
assume as well that random bonds may be antiferromagnetic, thus making connection to 
spin glasses. The discussion of these subjects, however, is out of the scope of the present 
exposition. 
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4 The N-replicas (lattice) 



Within the iV-replica scheme, see (f|), we have to average the replicated partition function. 
From ([TDD, the replicated partition function is: 

N N N 

Q N {y}= / II X{d-C n dC n X{[5{l-y mn )X{e s -+5{y mn ) ![«„], (12) 

a=l mn mn a =l a=l 

where S£ n (a = 1,...,N) are the replicas of the RS-pure local action from (|TT|). The 
averaging over the disorder with distribution (Eh then yields: 



N N N 



r " - a 
mn mn 



Q N {y} = I II UdCmndC n u[pu e 5 ™ + a -p) n 

a=l mn mn a=l a=l 

r N N 1 — p N 

/ II II dc mn dc mn exp ^2 [ 'S'mn "I IT C mn^mn eX P ( — ^mn) 

a=l mn mn a=l P a=l 

/ IT IT dc mn dc mn exp ^2 [ ^mra "I ~ IT C mn^mn eX P (^mn) ■ (1^) 



r2 



r 2 



a=l mn mn a=l P a=l 



We elaborate the product from the first line into a unique exponential making use of the 
nilpotent property of fermions, (c^„) 2 = (c^„) 2 = 0. In the final line, referring to the same 
property, we replace -S^ n by the reduced form, A^ n = i^^c°_ ln c° n _ 1 , taking into ac- 
count that prefactor c^ n c^ n annihilates the correspondent fermions in exp (— S^ n ). Yet an- 
other possibility, preferably to be used below, is the choice: A™ n = t^tff^dmC^n) (d n Cmr 



where d m , d n are lattice derivatives: d m c« n = c« n - c«_ ln , d n c« n = C - 



We have thus formulated the exact lattice fermionic theory for the 2D Ising model with 



quenched site dilution. The integral resulting in (|T^) can be written shortly in the form: 

Q N {y} = I D exp (So + A S int ) , A = ^— ^ , (14) 
J v 



where D is the replicated fermionic measure, So is the Gaussian part of the action, which 
is the replicated action for the RS pure case, eq. (pT|) , while the non-Gaussian term of 
interaction can be written, in particular, as follows: 



Ao 'S'int Ao ^2 { IT C mn C mn eX P (^mn^rnn i^m C mn) (^n C mn) )} 
mn a=l 
N 

= Ao ^2 { IT C mn^mn 1 + ^mn^mn (^m c m n) (^n^mn) } ■ (1^) 



mn a=l 



The interaction appears to be of order 2N in fermions, where N is the number of replicas, 
with N — > at final stages. This is a somewhat unusual feature of the resulting model. In 
a sense, in ([15]), at each lattice site we have the product of 'holes' (virtual dilute sites) over 
the replicas. Respectively, the coupling constant is of order 1 — p, which is a fraction of the 
dilute sites over the lattice. Despite some unexpected features, the interaction ( |T5|) is of a 
relatively simple and compact form. The resulting fermionic theory can be used, probably, 
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in future research, for the analysis of strong and moderated dilution (it might be important 
to take properly into account short lattice distances). For weak dilution near T c , where the 
continuum- limit approximation can be applied, the standard four-fermion term in fact arises 
in the effective action due to the interplay of higher- and low-momentum lattice fermionic 
modes. This is discussed in section 7. Before, we comment shortly on the lattice analytics 
and continuum limit for the pure 2D Ising model. In what follows we assume ferromagnetic 
case with the homogeneous distribution of bond coupling parameters over the lattice. 



5 The pure case (lattice) 

From flnp specified to the homogeneous distribution of bond couplings, we obtain the Gaus- 
sian integral for the partition function of the regular 2D Ising model: 



ran ran 

*1*2 C m — lijCmn- 1 t\ C mn C m —\ n -\- t 2 ^mn—X^ran ] > (^^) 

where t\ j2 = tanh bi )2 , bx t2 = J\,i/kT. Let us pass to the momentum space by Fourier 
substitution: 



u pq c i ^mn nr~? pq ' v / 



•'ran f T 'I ' J **** ' ma / r 2 

V J-> pq V J-' pq 



with —L/2 < p, q < L/2, where L is lattice length. The partition function becomes: 



Q= I Y[dc pq dc pq exp [(l-he^ -t 2 e i2 t" -t^e' 



2ttj] i ■ 2-n-g . 

L L ) c„„c r 



pq pq 

pq pq 

■ 2irp ■ 2irq 

+ t 1 e l — c L _ pL _ q c pq + t 2 e l — c pq c L _ pL _ q 

Making the pq <-» L—pL—q symmetrization of the sum (needed since variables with momenta 
pq and L — pL — q interact) we find the integral (p~8[ ) to be decoupled into a product of simple 
low-dimensional integrals, each over dc pq dc pq dcL_ P L-qdcL- P L-q- The pg-integral factor can 
be readily evaluated making use of the elementary rules, see ([]), which results in the explicit 
solution for Q 2 . For similar momentum-space calculations see also |29j. Taking then the 
limit of infinite lattice, for the free energy per site we obtain: 



7T 7T 



-PfQ = \l I ln ft 1 + + tl) - 2ti(l - 4) cosp - 2*3(1 - t?) cosgl , (19) 



2 J J (2vr) : 

which is equivalent to the Onsager solution As it follows from the analysis of the free- 
energy form, in the ferromagnetic case, t±, t 2 > 0, the critical point is given by the condition: 

1 - h - t 2 - t x t 2 = . (20) 

The specific heat exhibits the logarithmic singularity near the critical point, C ~ A c \ In |r||, 
t ~ |T — T c | — > 0, where A c is the critical amplitude. For the symmetric lattice (ti = t 2 ) the 
amplitude is: A c = (8/V) where b c =\ ln(l + a/2) = (J/kT c ) is the dimensionless inverse 
critical temperature. 
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6 Majorana and Dirac Fields 



Let us make a comment on the continuum-limit interpretation of the pure 2D Ising model 



near the critical point, also see [30, 31]. It is interesting that already the exact lattice action 
from ([16|) can be written in the field-theoretical like form, cf. ||31|| . In terms of the lattice 
derivatives d m c mn = c mn - c m _ ln , 8 n c mn = c mn - c mn _i, this action becomes: 

S (c, c) p Ur e ^ ] |^ Tfl Cmn Cmn ~\~ Al (<9 m C mn )c mn A 2 C mn (8 n C mn ^ 
mn 

~\~ t\ {dmCmn)C mn ~\~ t 2 C mn (8 n C mn } ^1^2 (^mCmn) (^8 n C mn j , (21) 

with m = 1 — ti — t2 — t\t2, \i = t\ + t\t2 , A 2 = t2 + t\t2- Evidently, m plays the 
role of mass, the critical point corresponds to m = 0. Taking the formal limit to the 
continuum space: mn — > (xi,%2) = %, d m ,d n — > 81,82, and changing the notion for the 
fields: c mn , c mn — > 4> (%) , 4> (x) , we obtain the continuum-limit counterpart of the exact lattice 
action (^l|) of the same formal structure. The second order kinetic term with 8182 can be 
ignored. The remaining action is of the Majorana type and can be transformed into a 



canonical form by a suitable rotation and rescaling of the fields, if),ift — > ipi,ip2, cf. |30|, plf . 
In this way we obtain the action in a canonical Majorana form: 

Styi,ik) = J d 2 x[m^2+ ^iU d i + id ^i+^U- d i+ id M, (22) 
with the rescaled mass: 

l — ti—t 2 — tit2 , . 

m = , (23) 

V 2 (tit 2 ) c 

where ( ) c stands for the criticality condition, eq. (p0|). 

The 2D Majorana action (|22]) can be rewritten also in the matrix notation. Assuming 
°2) °"3 to be the standard Pauli matrices, we find: 

S maior = l Jd 2 x^(x)[m+8]^(x), ^(x) = ^ T (x)(*cr 2 ), 9 = 7^+72$, (24) 

where ^ and \I/ are the two-component left and right Majorana spinors, \I/ = (ip 1,1^2), an d 
7i — a \i 72 — < J 2 are the 2D gamma-matrices. Finally, we can pass to the Dirac theory 
of charged fermions by doubling the number of fermions in the Majorana representation 
14] , |31 |. Taking two identical copies S' and S" of the Majorana action, eq. (24), we consider 



the combined action £ dirac = (S' + £") major- Introducing new Dirac fields by substitution 
= -j=(ty ), = 4g — ity ),we find the combined action in the form: 

5 dirac = J d 2 xy{x)[m +8]V (x) , (25) 



where m and 8 are the same as in (pip, while and \I/ are now charged Dirac spinors with 
four independent components. 
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The formal meaning of the transformation to the continuum limit can be most readily 
understood in the momentum space. The change d m ,d n —>■ di,&2 in the momentum space 
corresponds to the approximation: 

e w — i ~ ip ( e «<? _ i ~ j,q , (26) 

where p, q are the quasidiscrete Fourier momenta, p,q ^j 2 , — 7r < p, q < tt, see (0), 
and |p| below is modulus of (p,q). The approximation like (p6|) is legal only in the neigh- 
borhood of the origin in the momentum space, < \p\ < ko, where ko -C 1 is some fixed 
small cut-off momentum. There is also an effective infrared bound for \p\ introduced by the 
mass, which implies that m is small and we are near T c . The continuum-limit formulation 
thus captures the essential features of the exact lattice theory in the low-momentum sector 
responsible for the singularities in the thermodynamic functions and the large-distance be- 
haviour of the correlations as T — > T c (m — >• 0). In this region, the approximation is exact. 
In particular, the logarithmic singularity in the specific heat can still be recovered from ( ^4|) 
or (|25|), with the exact value of the critical amplitude. 



7 Gross-Neveu model (N —> 0) 

To extract the continuum-limit theory for the model with interaction dl3|)-(|T5[), let us distin- 
guish explicitly the lattice modes, cf. (|T?D, with higher and low fermionic momenta, \p\ > ko 
and |p| < k , where /c 1 is small fixed cut-off. Introducing separable notation a, if) for the 
corespondent fields, in the coordinate space we write: 



The Gaussian part of the exact lattice action (|14D is additive in the momentum space, so 
So = S a + Sip, where S a and S^p are independent. However, the non-Gaussian term of 
interaction (|i~5|) is non-additive with respect to c — ► a + if). We have thus to substitute 
explicitly c — > a + if) into this term, and then integrate over the higher-momentum lattice 
modes (a mn ) \ P \>k - This is done below for weak dilution, in the first order of perturbation in 
A ~l-p<l. 

Writing the action ([14]) in the form S = S a + S^ + \o Si nt (a, ip), An = — -, in the first 
o 



V 

order in An we obtain: 



Q N ~ J DaD^e^+^il + XoSint (a,ifj)) 



Substituting fflj ) into (Si nt (a, if>)) a , we follow the standard field-theoretical prescriptions for 
critical fluctuations that only the most important quartic and quadratic terms in the if>,ip 
fields are actually to be taken into account [14], [15| . In terms of the characteristic fermionic 



8 



averages (A) and (B) presenting the contribution of the higher-momenta lattice modes, we 
find: 



N 



S int (a,i>)) a = EII [(A) + (B) + - 

van a=l 



0=1 



(29) 



where 



(A) 
(B) 



{aa exp (tit 2 did d 2 a)) a = (aa) a + t\t 2 (act d\a 82a)) a 
(exp (ht 2 dia d 2 a)) a = 1 + ht 2 (diad 2 a)) a . 



(30) 



These are the lattice averages for pure model (16) with infrared truncation at k . Here we 



assume the evident abbreviation with dropped ( ^ nn index and 8%, d 2 «-> d m , d n correspondence. 
For instance, (diad 2 a) a means (<9 m a^ n <9 n a,^ n ) a , etc. The dependence on small momentum 
cut-off k can be in fact ignored (put k = 0). The averages are then to be evaluated with 
the pure lattice action, eqs. (|16] ) - ([18|) , and are ^ n independent. 

In (|29|), in the first line (...) is the abbreviation for the 'gaussian kinetic terms', which can 
be ignored, and in the second line (...) stands for the total sum of all irrelevant corrections, 
including the gaussian kinetic terms, quartic kinetic terms, and higher-order corrections 
in the 0>, $ fields of any kind. The effective action for the ■0,0 fields is S c s = (SVOfree + 
A (S int (a, if))) a . The terms preserved explicitly in ( ^9]) contribute to the Gaussian part 
of this action (modification of mass) and generate the four-fermion interaction. A trivial 
constant term may be ignored. The Gaussian part of the action S c q appears in the form 
of the replicated continuum-limit version of the action fl2T|), cf. (|22|), but with a modified 
mass. After uniformization of the Gaussian part and the correspondent rescaling of the fields, 
i/jj't/j — > 0i, i/) 2 , we obtain the effective continuum- limit action in the form of the iV-colored 
Gross-Neveu model || P| : 

5 G _ N = J d 2 x{ E [m N ^2 + 2 W & + W 



a=l 



4* 



N 2 

di+id 2 )^] +^[E WW] }, 



with the effective mass and coupling constant given as follows: 

, /A , N l-p (B) 



m N 



9n = (A) 



N 



2{tit 2 ) c 
l-p(B) 2 1 



(A) 



P (A)M(M 2 )c 



(31) 



(32) 



(33) 



Since the replica limit N — > is assumed at final stages, one can put N = and (A) 
already in (|3~2"D and (0). In the matrix notation, the action (RTT) becomes: 



N 



1 N r r ~\ N 

-E /rf 2 x^(a:)(m + 9)vl>-(x)+^ d 2 x [- E * a l 



5 



G-N 



x)^ a (x) 



a=l 



(34) 
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The Gaussian part here is the replicated Majorana action (2jj) with modified mass (|32|). 
By analogy with the pure case, we can pass as well to the Dirac field interpretation. The 
doubling of Majorana fermions can be realized simply by doubling the number of replicas, 
N -> 2N, and we find: 



S 



N 

<;-N = 

a=l 



d 2 x^ a (x)(m + d)^ a (x)+g d 2 x 



N 



(x)V a (x) 



a=l 



(35) 



where ty a ^ a are charged Dirac spinors, cf. (P5|). In (|34|), (j35|) we assume the mass and 
charge already at N = 0. 

The N = Gross-Neveu model was first introduced as an effective field theory for weak 
bond dilution in the pioneering papers by Dotsenko and Dotsenko (DD) Their theory 

has been reanalyzed later on as regards to the behaviour of correlations by Shalaev, Shankar, 
and Ludwig (SSL) 0-0, see WM for a comprehensive review. An advanced analysis of 



correlations in the N = G-N model in terms of the perturbed conformal field theories has 
been performed quite recently by Dotsenko, Picco and Pujol [[JTJ, |34|j. We thus expect the 
same changes in critical behaviour inspired by weak site dilution as those predicted in the 
DD-SSL theories for weak bond dilution. From the renormalization-group (RG) analysis of 
the N = G-N model it follows that in a closed vicinity of the critical point the singularity 
in the specific heat is doubled-logarithmic || |15| : 



C oc — In f 1 



^,1 

7T T 



(36) 



| that the correlation radius gains the logarithmic correction, 
~ L in the pure case, where r ~ \T — T c \). The two-spin 



The RG calculation also shows 
R c ~ r-^lnr] 1 / 2 (cf. R c ~ r 
correlation function G (R) = (o~(Q)o~(R)) taken exactly at T c remains to be the same as for 

-1/4 



the pure system |14]: G (R) ~ R' 1 ' 4 , R — > oo. [That no logarithmic or other corrections 
appear due to the disorder in G(R) in the N = G-N model can now be assumed to be 
definitely established |14|, |33| , p4| , this question was a subject of some debate in previous 
decade 0-0]- Assuming the standard scaling arguments to be true, this then implies 
that the magnetization and magnetic susceptibility gain the logarithmic corrections to the 
pure-case power singularities at criticality [|TJ, E2J: M (r | 0) ~ |t|^ 8 | lnr| _1//16 , x (r | 0) ~ 
|r| _7 / 4 | lnr| 7 / 8 , M(0 | h) ~ h 1 ^ 15 , where r ~ |T — T c \ and h is a small magnetic field. These 
predictions have been checked also, for weak dilution, in a series of the precise Monte-Carlo 
experiments fl7|l , fL8| , pT| , [Q. The question about the effects produced by strong and 
moderated dilution is less evident, however, both on theoretical [IT 7 !, Q and experimental 



(Monte-Carlo) |20], 55] sides, this holds for both bond and site dilution. 

In conclusion to this section, let us specify in more detail the parameters of the N = G- 
N model for site dilution for the isotropic lattice (ti — t 2 , or J\ — J 2 ). This will yield also the 
exact value for the initial slope of the T c (p) curve. The effective mass and coupling constant 
include the lattice constants (^4) and (B). Since the dilution is weak, and we are near T c of 
the pure system, it is reasonable to specify these parameters exactly at T c of the pure system. 
Also, the dependence on small cut-off momenta, ko <C 1, can be ignored. The averages (A) 
and (B) are then to be taken for the pure lattice system at T c . If t\ ^ t 2 , however, these 
parameters still remain be functions of lattice anisotropy. They become definite numbers 
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in the isotropic case. By a straightforward though somewhat lengthy calculation we have 
obtained the following values for the isotropic lattice (J\ — J2, T — T c , k = 0): 



(B) c = 2t c (l + —= --) = 1. 150517, t c = V2 



V2 



71 



^r 1 = 4 V2t c (l + — ) 1 = 1.615786. (37) 
{A) c V TV J 

For the isotropic lattice, the critical point is given by condition (1 — It — t 2 ) c = 0, with the 
solution t c = e~ 2bc = V2 — 1, b c = \ ln(l + a/2). The pure-case mass m = m(l) ( ^3|) then 
appears in the form: m(l) ~ 2{t c — t)/t c ~ 46 C (T — T c )/T c . Taking into account (|3~7D , the 
effective mass (N = 0) of the G-N model becomes: 

T — T 1 — v 4 
m = 4 6 C - c + ^— (38) 

= (1.763) ^-^ + (2.758) 

c P 

while the coupling constant (N = 0) appears in the form: 

= — 7VL 8 /o/ ^2 - ( 3 - 804 ) — ' ( 39 ) 

P (1 + V2/7T) 2 P 

These are the values of mass and charge of the iV = Gross- Neveu model (|31"D-(|3"5D for weak 
site dilution near T c in the isotropic case. These values are exact in the linear orders in 1 — p 
and T — T c . 



The expression for mass fl38|) enables one to define the critical slope for the T c (p) — p curve 
at p = 1. The initial slope is the coefficient in the expansion T c (p) = T c (l) [ 1 — S c — - + ...]. 
Dividing both sides of the equation m (p) = by 1 — p and taking the limit p — > 1, we find: 



T c dp 



p=1 1 + V^/tt 



-1 

1.564785, (40) 



which is in agreement with the earlier calculation of S c [P5(| . As distinct from a direct spin- 
lattice perturbative analysis in R3] , the present calculation includes the N- replica trick with 
fermions. 



8 Concluding discussion 

We have applied a new noncombinatorial method of fermionization to clarify the fermionic 
structure of the 2D Ising model with quenched site dilution. In fact, the method works 
equally well for both site and bond dilution. We have specified the discussion merely to 
the case of site dilution as being less studied theoretically. For the first step, the partition 
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function of the 2DIM with fixed site dilution and arbitrary inhomogeneous distribution of 
the bond coupling parameters over the lattice was transformed into a fermionic Gaussian 
integral, eq. @. The resulting integral was then simplified integrating out extra fermionic 
degrees of freedom. This yields the Gaussian like representation with two fermionic vari- 
ables per site for fixed site dilution, eqs. (|9|)-(|TT1). Even for the pure case, this reduction 
of extra fermionic degrees of freedom essentially simplifies the analytics of any kind and 
also illuminates Majorana-Dirac structures in the pure 2DIM already at lattice level. The 
quenched averaging over the disorder in the iV-replica scheme results, in turn, in the ex- 
act lattice theory of interacting fermions for the case of site dilution, eqs. fli~^ ) - (|i~5|) . The 
interaction in this exact lattice theory appears to be of order 2N in fermions, where N is 
the number of replicas. However, the continuum-limit approximation for weak dilution near 
criticality results in the standard N = Gross-Neveu model with four-fermion interaction. 
The N = G-N model has been already intensively analyzed as an effective field theory for 
bond dilution || |TJJ . This then implies similar behaviour of the thermodynamic functions 
and correlations in the disordered 2D Ising model both for site and bond dilution near T c , 
at least for small amount of impurities. The predictions of the N = G-N model for weak 
dilution are the log-log singularity in the specific heat and log- correct ions in M and \ , see 
section 7 for more comments. We have also specified in detail the parameters (mass and 
charge) of the N = G-N model for site dilution in terms of the lattice fermion Green's 
functions and evaluated the initial slope of the T c — p curve. The present analysis thus 
supports the hypothesis about the universality in critical behaviour produced by different 
sorts of disorder for a small fraction of impurities in the ferromagnetic 2D Ising model. 
The question about the universality between weak and strong dilution remains to be open, 
in essence. A more straightforward analysis, directly in lattice interpretation, seems to be 
needed. A compact form of interaction in the exact lattice theory (ITS)- (15) provides grounds 
to suppose further progress in this respect. Finally, let us note that the Gaussian fermionic 
representations with fixed disorder like (|8|), (^), ([□]) yet preserve all the information about 
the original spin- variable model ([3|) that is included in the hamiltonian (JJ). There are no 
also any restrictions, at this stage, on the sign of the bond coupling parameters. Therefore, 
these representations can be used as the starting point to try other problems with disorder 
in the 2D Ising model, beyond site or bond dilution (spin glasses, models with annealed 
or regular inhomogeneities, etc). As well, one can try to apply other methods of quenched 
averaging, beyond the iV-replica scheme. The later possibility may be favorable concerning 
the problem of strong and moderated dilution. 
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